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(— I ■ Abstract: We investigate continuous Hawking-Page transitions in Einstein's grav- 

ity coupled to a scalar field with an arbitrary potential in the weak gravity limit. 
We show that this is only possible in a singular limit where the black-hole horizon 
Q ■ marginally traps a curvature singularity. Depending on the subleading terms in the 

^ ■ potential, a rich variety of continuous phase transitions arise. Our examples include 

O ' second and higher order, including the Berezinskii-Kosterlitz- Thouless type. In the 

t^^ ■ case when the scalar is dilaton, the condition for continuous phase transitions lead 

Q ■ to (asymptotically) linear-dilaton background. We obtain the scaling laws of ther- 

modynamic functions, as well as the viscosity coefficients near the transition. In the 
limit of weak gravitational interactions, the bulk viscosity asymptotes to a universal 
^ ■ constant, independent of the details of the scalar potential. As a byproduct of our 

c^ ' analysis we obtain a one-parameter family of kink solutions in arbitrary dimension 

d that interpolate between AdS near the boundary and linear-dilaton background 
in the deep interior. The continuous Hawking-Page transitions found here serve as 
holographic models for normal-to superfluid transitions. 
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thermodynamics. 
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1. Introduction 

Thermodynamics of radiation in asymptotically Anti-de- Sitter spaces is quite inter- 
esting. It was shown by Hawking and Page in |jl| that, as one heats up a gas of 
gravitons in the global AdS space-time, one finds a phase transition into nucleation 
of black-holes at a finite temperature T^ ~ i~^ where i is the AdS radius. At low 
temperatures, the only equilibrium solution to Einstein's equations with asymptot- 
ically AdS boundary condition is thermal radiation. Above a certain temperature 
Tmin ~ ^~^ a "small black- hole" and a "big black-hole" saddle coexist with the ther- 
mal gas (TG). Right above Tmin these solutions have free energy larger than the 
radiation, therefore they eventually tunnel into thermal radiation. However, above 
a certain value Tc > Tmin the big black-hole solution becomes the global minimum 
and the thermal gas instead becomes quantum mechanically unstable against tun- 
neling into the big black-hole.^ This is the so-called Hawking-Page transition in 
asymptotically AdS space-times. 

Witten in proposed that the Hawking-Page transition in the AdS^ x S^ solu- 
tion of IIB string theory corresponds to a confinement-deconfinement type transition 
in J\f = 4 Yang-Mills theory on S^ with radius R ^ i, through the AdS/CFT 
corespondence0]. There is a transition at some Tc ~ R~^ ~ i~^ and it only occurs 
in the large color limit N -^ 00, where the number of degrees of freedom becomes 
infinite. A more interesting case is when the field theory is defined on R^, which can 
be obtained as the large radius limit of the sphere i? — ?■ 00. Clearly Tc — ?■ 0"^ in this 
case. The corresponding statement in the dual field theory is that the A/" = 4 sYM 
theory is conformal, and when placed on S^ x S^ with S^ being the time-circle with 



^Thc small BH solution has negative specific heat and it is always subdominant with respect to 
the big BH and TG saddles. It will play no role in what follows. 
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radius 1/T, the thermodynainic quantities can only depend on the dimensionless ra- 
tio TR. Therefore, there cannot be any interesting thermodynamics in the flat space 
limit. 

The situation changes, when there is an additional intrinsic scale in the problem 
that does happen in non-conformal theories. A basic example is pure Yang-Mills 
theory with gauge group SU{N). The theory has running coupling, hence a dynam- 
ical scale A is generated through dimensional transmutation. Indeed, lattice studies 
provide strong evidence for a second order deconfinement transition for N = 2 and 
a first order deconfinement transition for A^ > 2 in 3 space dimensions roughly at 

By a natural generalization of the ideas in 0, one expects that Hawking-Page 
type transitions in asymptotically AdS space-times dual to Yang-Mills type non- 
conformal theories^ should correspond to confinement- deconfinement transitions in 
these theories. In , it was shown that the simple set-up of Einstein gravity coupled 
to a scalar field $ with some potential V{^) allows for asymptotically AdS type 
solutions with a non-trivial intrinsic scale A that is generated due to the running 
scalar field. It was proposed in P that the ground state (zero T) solution generates 
color confinement in the dual field theory, if and only if the IR asymptotics of the 
potential obeys, 

y(<l>)^e3*$^, P>0, as $^00, (1.1) 

in 5 dimensions. A background that solves a potential of the form (|1.1|) in the IR and 
becomes asymptotically AdS in the UV indeed exhibits a Hawking-Page transition 
that corresponds to the deconfinement transition in the dual gauge theory. 

In fact, the thermodynamics in the generic case P > is very similar to the 
global AdS case that we described above: There is always a thermal gas solution 
saddle at every T. This is the only equilibrium configuration up to a certain T^m- 
Above Tmin there are two black-hole solutions in addition to the thermal gas: the 
small and the big BHs, just as in global AdS. At a certain Tc > Tmm the big BH 
solution becomes the global minimum of the theory. This temperature corresponds to 
the deconfinement temperature in the proposed dual YM theory, where the glueballs 
deconfine their glue contents. For P > 0, one can quite generally show that, there 
is a finite latent heat and the transition is of first order. 

In this paper, we ask the general question, "under what conditions gravity ex- 
hibits continuous Hawking-Page type phase transitions?" By the word continuous, 
we mean second or higher order transitions where the entropy is continuous across 
the transition. In other words, the latent heat vanishes. Regardless of the matter 
content of the gravity theory, the basic condition for a continuous HP transition can 
easily be obtained in the classical limit Gd — ?■ 0. In this limit we arrive at the (more 



^By a Yang-Mills type theory we mean a gauge-theory that exhibits quark confinement in the 
ground state. 
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or less obvious) fact that, a continuous transition in this limit can only occur in a 
singular region of space-time, where the BH horizon marginally traps a curvature 
singularity. This is easily understood: in the weak curvature limit, the contribution 
of thermal radiation to the entropy is negligible, thus the entropy of the TG solution 
is very small (zero in the strict Gd — ^ limit). Then, continuity of entropy requires 
that the horizon size of the BH solution becomes very small as one approaches the 
transition temperature from above. 

To be more definite, we approach the problem of continuous transitions in the 
context of Einstein's gravity coupled to a scalar. In fact, it was already observed 
in [^ that, the aforementioned HP transition may be continuous in the marginal 
limit P = in ( |1.1D . Here, we focus in this particular case in more detail and also 
generalize the analysis to arbitrary d+1 dimensions. Our general result can be stated 
as follows: there exists a continuous Hawking-Page transition in d + 1 dimensional 
Einstein-scalar theory if and only if 
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^($) ^ e—"" (1 + Vsu^m , as $ ^ oo, (1.2) 

where Vsub denotes subleading corrections that vanish in the deep IR. We also find 
that these corrections determine the type of the transition: For an exponential fall- 
off in Vsub, one finds nth order transitions for n > 2. In this case the free energy 
difference scale as, 

AF ~ (T - Te)", as T^ T+ . (1.3) 

On the other hand, for a power-law fall-off, a Berezinskii-Kosterlitz-Thouless type 
scaling arises: 

AF~e-^*"", asT^T+, (1.4) 



where a > 1 is a constant. 

A good indication that the marginal case P = 0, corresponds to continuous tran- 
sitions is obtained by investigating the equilibrium solutions of Vconf- As mentioned 
above, for P > generically there are three solutions at high temperatures: the ther- 
mal gas, the big and the small black-holes. In the marginal limit P — ?► the small 
black-hole solution disappears. Therefore, in the high T regime there exist only the 
TG and the big BH branches. In fact, in this case Tmin and T^ can be shown to coin- 
cide and the thermodynamics becomes much simpler than the generic P > case: 
The low T (high T) phase is dominated by the TG (big BH) saddle and these phases 
are separated by a Hawking-Page transition at T^. This is precisely what one expects 
in a continuous phase transition: At T^ the new saddle is created and immediately 
becomes the global minimum of the free energy. Indeed the small BH solution can be 
considered as a kink that connects the TG and the big BH solutions which can only 
exist in first order transitions. In the limit where the transition becomes continuous 
this kink should disappear. 
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As we mentioned above, in the classical limit, the continuous transitions gener- 
ically occur in a singular region where the horizon area becomes very small. In 
the particular example of the Einstein-scalar gravity with potential ( p. .11) and with 
P = 0, this region corresponds to $/i — )• oo where $/i is the value of the scalar on the 
horizon^. In this limit, it was shown in that the BH and the TG solutions become 
degenerate, that is indeed another basic property of second or higher order transi- 
tions. In fact, in the strict classical limit Go — )■ 0, the transition point corresponds 
to $£ = oo. 

Clearly we cannot trust Einstein's gravity in this singular region. On the other 
hand we would like to study the scaling of thermodynamic functions with T — T^., 
an information contained in this singular region. We avoid this problem here by 
fine-tuning the sub-leading term in (|1.2| ) such that the region T <Tc corresponds to 
a wide range $/i G ($s, oo) where $« is the value of $ where the scaling behavior in 
( |1.3| ) or ( |1.4| ) sets in. This is clearly a crude method that we can resort to, only in a 
phenomenological approach where we use the freedom in the choice of the potential to 
hide the problem of singularity. This method serves our purposes in this paper, that 
is to work out the qualitative behavior of thermodynamic functions in a continuous 
HP transition. However, a more detailed future study should take into account the 
a' corrections by investigating continuous HP transitions in the full string theory. 

Studying the IR geometry of the solutions to ( |1.2| ) we arrive at another interesting 
result: the BH and the TG geometries both asymptote to a linear-dilaton geometry, 
in the case when the scalar field is the dilaton. More precisely, the ground state 
(TG solution) that uniquely follows from the potential ( |1.2| ) is asymptotically linear- 
dilaton in the deep interior of the geometry. At the same time, the BH geometry 
approaches to the TG geometry in the limit T — )■ T^. Therefore both of the saddles 
are of linear-dilaton type in this region of space. It is well-known that the linear- 
dilaton geometry is an a' exact solution to the full string theory. This observation 
may prove important in studying the analogous problem in the full string theory. 

In the limit where we can ignore gravitational interactions and higher derivative 
corrections, Einstein's gravity predicts certain universal results around a continuous 
transition. One such quantity is the value of the bulk viscosity per entropy as T — )■ Tc. 
We find that this value is constant across the transition ^: 

-^ . J ,y asT^T,. (1.5) 

s lTi[a — 1) 

The result is universal in the sense that it is independent of the details of the theory 
such as the choice for the scalar potential ^($), the transition temperature T^ etc., 
and it is only subject to a' corrections in a gravity theory with higher derivative 
terms. 

^It can also be viewed as the horizon location in a coordinate system where the radial direction 
is given by $. 

^See S for another example where C,/ s is found to be constant. 
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We should mention that ahhough this work can be regarded as an investigation 
purely in gravity, we have holographic applications in mind. In this paper, we shall 
not specify the expected dual field theory. In fact most of our findings will not 
depend on the UV geometry either. Yet, we would like to have an asymptotically 
AdS geometry in the UV, to be able use the prescriptions of the AdS/CFT duality in 
holographic applications. A particularly interesting such application is studied in a 
companion paper ||^, where we propose that the continuous phase transitions found 
here are holographically related to normal-to-superfiuid transitions in spin-models. 

In the course of searching for backgrounds with AdS asymptotics, we also obtain 
a one-parameter family of analytic kink solutions that fiow from an AdS extremum 
to the linear-dilaton background in the string frame. 

In the next section, we list the general conditions for the Einstein-scalar theory 
to exhibit a continuous Hawking-Page transition. In section 3, we present the classifi- 
cation of the transitions that follow from different fall-off behavior of Vgub in ( [I.2[) . In 
this section, we also present a coordinate invariant method to study thermodynamics 
in Einstein-scalar gravity, that was first constructed in [0. Here we generalize this 
method here to arbitrary dimensions. Section 4 is devoted to derivation of the scal- 
ing behavior of thermodynamic functions: the free energy, the entropy, the specific 
heat and the speed of sound. In section 5 we study the shear and bulk deformations 
around the black-hole solution and study the viscosity-entropy ratios rj/s and (/s 
in the transition region. In section 6 we construct the aforementioned analytic kink 
soutions. Finally in section 7, we briefiy mention what type of corrections are ex- 
pected in our findings when the set-up is embedded in full string theory. Finally we 
present a discussion on our results and an outlook for future directions in the last 
section. Six appendices detail our calculations. 

2. Weak phase transitions and small black- holes 

2.1 Generalities 

In this section we determine the requirements for presence of continuous phase tran- 
sitions in the Einstein-scalar theory in c? + 1 dimensions, defined by the action: 

A = ^^ I df'+W^ [R - ad^^Y + V{<^)) + G.H. (2.1) 

We keep the normalization of the scalar kinetic term ^ > as a free parameter 
for later convenience. The last term, that we shall not need to specify here, is the 
Gibbons-Hawking term on the boundary. Go is the Newton's constant. The total 
number of bulk dimensions is denoted hj D = d + 1. 

There are only two types of backgrounds at finite T (Euclidean time is compacti- 
fied), with Poincare symmetries in d — 1 spatial dimensions, and and additional U{1) 
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symmetry in the Euclidean time direction: the thermal graviton gas, 

ds^ = e2^«('^) [dr^ + dxl_^ + de) , $ = <l'o(r), (2.2) 

and the black-hole, 

ds" = e'^^'-) {r\r)dr^ + dxl_^ + rftV(^)) , ^ = ^{r). (2.3) 

We define the coordinate system such that the boundary is located at r = 0. For 
the potentials V that we consider in this paper, there is a curvature singularity in 
the deep interior, at r = r^. In ( p.2|) , r runs up to singularity r^. In (|2.3| ) there is a 
horizon that cloaks this singularity at rh < Vg where /(r/j) = 0. t is the Euclidean 
time that is identified as t ~ t+l/T. This defines the temperature T of the associated 
thermodynamics. In the black-hole solution, the relation between the temperature 
and r/j is obtained in the standard way, by assuming that the near horizon geometry 
is conical and the conical singularity is absent: 

47rT = -/'(r,). (2.4) 

This identifies T and the surface gravity in the BH solution. 

The geometries ( p.2| ) and ( pl3| ) are the only backgrounds with the boundary con- 
dition that A{r) — )■ Ao(r), $(r) — )■ $o('") s'-iid /(r) — )■ 1 near the boundary^. We are 
interested in determining the necessary and sufficient conditions for existence of a 
continuous Hawking-Page type transition between (|2.2|) and ( p.3|) . We define a con- 
tinuous phase transition as the one where the entropy difference vanishes AS{Tc) = 
at the transition. This defines phase transitions of second and higher order, including 



infinite order, such as the BKT transition 10 



Thermodynamics is determined by the Gibbs free energy that is given by the 
on-shell value of the action { \i.l\ ). Taking the difference between the on- hell values 
on the BH and the TG solutions removes possible divergences due to infinite volume 
of the geometries and yields the regularized action: 

T-^AG{T) = A{BH) - A{TG) + ■■■ (2.5) 

where the ellipsis refer to corrections that are 0{Go) suppressed with respect to the 
first term. 



^This is somewhat rough, see M for precise matching conditions in the variable r, near the 
boundary. The easiest way to implement the matching is to consider a coordinate frame where $ 
is the radial variable and identifies $ = $o in the entire range of the two solutions. Then the only 
nontrivial condition is A{^) -^ Ao{^) and /(<&) — > 1 near the boundary. This is what we will do in 
section |[ 

^See ||ll|] and ||l^ for recent papers where the BKT scaling was observed in the context of 
condensed matter applications of holography. These papers consider quantum phase transitions 
however, where the scaling is not in T but in a quantum coupling. 
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The Newton's constant in ( p.l|) is important in determining the overall size of 
the thermodynamic functions. It also plays an important role in the correspondence 
between gravity and a dual field theory, where Gd should be proportional to the 
number of degrees of freedom in the dual quantum field theory. This is most clear 
when the proposed dual theory is is a gauge theory with A^ colors: 

^ M^-^N\ (2.6) 



WtvGd 



where Mp is a "normalized" Planck scale, that is generally of the same order as the 
typical curvature of the background^. On the gravity side the factor N arises in a 
similar fashion to the AdS^ x S^ case. Here it is supposed to be proportional to the 
RR d+1 form flux on the geometry [^] . As explained in [^ the space- filling flux can 
be integrated out by replacing it with its equation of motion in the action, and it 
only contributes to the potential V{^). The potential V that is written in (|2.1|) is 
supposed to include this contribution. This procedure is completely analogous to the 
5-dimensional case discussed in and we refer the reader to this paper for details. 

Clearly, in order to be able to ignore interactions between bulk fields one needs 
to take the large N limit. This can easily be done when the scalar is dilaton: in an 
asymptotically AdS solution where the boundary value of the dilaton is constant $0, 
the large N limit is given by sending the dilaton $o — > — oo, and the flux N -^ oo 
such that e'^°N = const. 

In this work, we will mostly be interested in the scaling of the various thermo- 
dynamic quantities with T and not in the overall size of these quantities. For that 
purpose, it is more convenient to factor out the overall Gd dependence by dividing 
the free energy by Vd-iN"^ where Vd-i is the volume of the spatial part, and A^^ is 
defined by (|2.6|). We, thus define the normahzed free energy density, 



AF(r) = ^^ ^^_^^, {A{BH) - A{TG)) + ■ ■ ■ ^ F,{T) + 5F{T). (2.7) 

Here Fq is defined as the leading 0{N^) piece in the free energy density, and 5F refers 
to the 0{N~'^) corrections. The large limit N ^^ oo {Go -^ 0) defines a saddle point 
approximation from the gravitational point of view. The corrections that we denote 



by 6F in (2/7) arise, for finite but small Gd from the determinant of gravitational 



fluctuations around the saddle: 

6F{T) = Tr log Obh - Tr log Otg, (2.8) 

where O is the d'Alembertian operator of the fluctuations around the saddle geome- 
tries.^ 



''For example, in the improved holographic QCD models Q, it is given by Mp = (457r^)^3^~-'^ 
where i is the curvature radius of the asymptotically AdS geometry near the boundary. 

^In principle, for finite Gd, one has to make sure that this difference is free of divergences. By 



Although, we shall mostly be interested in the leading Fq term in (2/7), we note 
the following important consequence of the SF correction: Using the 1st law, the 
entropy (density) difference is given by, ^ 



where we normalize the entropy the same way as in (p.7|) . The Stg term above has 
no leading 0{N^) piece because the TG geometry has no horizon. Therefore, one 
has, 

AS = Sbh-^. (2.10) 

We shall ignore the l/N"^ corrections throughout the paper and drop the second term 
above, hence the entropy (density) difference of the dual theory is entirely determined 
by the black- hole entropy that is therefore given by the horizon size at Vh- 

(d-l)A(rfc) 

A^ = Sbh = ^^ ^, = 47rMp^-ie('^-i)^(^''^), (2.11) 

where we used the relation ( p.6| ) in the second equation. 

Therefore, we learn from ( p.ll| ) that — up to l/N"^ corrections that we neglect in 
this paper — a necessary condition for a continuous phase transition is that the black- 
hole becomes infinitesimal at T^. This is an infinitesimal hlack-hole^^ with vanishing 
classical horizon size.^^ In the backgrounds we consider there is a single curvature 
singularity at r = Tc. Precisely at Tc, the singularity is marginally trapped by the 
black-hole, and T^ corresponds to the temperature where the horizon approaches the 
singularity r^ -^ r^ ^^. 

There are further conditions for the existence of a phase transition as explained 
in the following. Before we present the derivation in the next subsection, let us 
state the final result: In the Einstein- scalar gravity, there exists a continuous phase 
transition into nucleation of black-holes if and only if the scalar potential V{^) has 
the asymptotic behavior, 

y($) ^ e^\/^*, $ ^ cx). (2.12) 



the boundary condition that BH asymptotes to TG near the boundary, the two spectra in the UV 
are guaranteed to be the same (this can be checked by a WKB analysis) , hence the difference should 
be convergent. 

^The backgrounds that we consider in this paper satisfy the required energy conditions and the 
laws of thermodynamics should apply. 

^°Not to be confused with the aforementioned "small black-hole" saddle which is a black-hole 
with a negative specific heat. In all of the examples we consider in this paper, we only consider 
black- holes with non- negative specific heat. 

^^Note that the black- hole at Tj, should be non-extremal, but with vanishing enthalpy. 

"'^^One may consider the possibility of string corrections at this point. We will discuss such possible 
corrections in the last section. 
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This result is of course subject to field redefinitions in the action (|2.1| ). Below, we 
present a derivation where we stick to the form of the action given in ( p.l|) . 

2.2 Conditions for presence of continuous transitions 

We are interested in a Hawking-Page type phase transition when the thermal gas of 
gravitons ( |2.2|) becomes unstable against nucleation of black holes ( |2.3| ) and we ask 
the question under what conditions this transition is higher than first order? The 
thermodynamic requirements are clear: 

i. 3 Tc such that < T^ < oo and AF(Tc) = 0, 

ii. AS{T,) = 0. 

Solution of these conditions for a generic scalar potential V in ( p.l| ) may seem rather 
arbitrary. This is not the case however. As we show below, the conditions above 
imply the specific form ( 2.12| ) for the asymptotic form of V. 



Before solving the conditions above, first we have to define the various thermody- 
namic functions in terms of geometric quantities. We solve the equations of motion 
and derive the thermodynamics in a coordinate independent approach in section ^72. 



Here we shall first work with the coordinate frame — that we call the r-frame — defined 
by the metrics ( p. 21) and (|2.3|) for simplicity. In the r-frame, one derives the following 



Einstein and scalar equations of motion from (p.l|) : 



A" - A'^ + -^^'^ = 0, (2.13) 

f" + id-l)A'f' = 0, (2.14) 

{d - l)A'^f + A'f + A"f - T^e^^ = 0. (2. 15) 



One easily solves (p.l4|) to obtain the "blackness function" /(r) in terms of the scale 
factor as, 

/(^) = l- j;.e-(^-i)A - (2-16) 

Then the temperature of the BH is given by equation (|2.4]) : 

ri'h 

T~^ = ^T,e^d-l)A(r,) ^-(rf-l)AM^^_ (2.17) 

Jo 

The difference between the entropy densities of the BH and the TG solutions is given 
by the BH entropy density ( [^.IID up to l/N"^ corrections that we ignore from now 



on: 



AC — (d-l)A{Th) /o io\ 
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The difference in the free energy densities ( p.7|) can equally well be evaluated by 
integrating the first law of thermodynamics ( p.9|) , [0]:^'^ 



Here Vc is the value of the horizon size that corresponds to the phase transition, 
T{rc) = Tc, where the difference in free energies should vanish. We note that, depen- 
dence on r/i in the integrand of ( p.l9|) is two-fold: There is the explicit dependence on 
the variable Vh and there is the implicit dependence that arise because the function 
A{r) itself changes as Vh is varied^'*. As we showed in 0], this implicit dependence 
on r/j is suppressed in the limit r^ close to the singularity r^. This is a fact that will 
prove important in the analysis of the next section. 

It was demonstrated in that, for any domain- wall type geometry of the form 
( p.2| ) with the condition that the scale factor obeys A'{ro) < near the boundary ^^, 
the possible asymtotics in the deep interior can be classified as: 

A. An asymptotically AdS geometry in the deep interior as r — )■ cxd, 

B. A curvature singularity at finite r = r^ where (J)q ^ oo and Aq — t- — oo, and 

C. A curvature singularity at r = oo where 0o — ^ oo and Aq — )• — oo. 

We want to determine which type of geometry above admits a Hawking-Page tran- 
sition by satisfying the conditions i and ii stated in the beginning of this section. 

We first focus on the condition ii: For a continuous phase transition, one should 
require that the entropy difference ( |2.18| ) should vanish at the transition point r^- 
Then, from the classification above we learn that r^ should correspond to either of 
the singular points r = r^ in case ii, or r = oo in the cases i or iii, because that is 
the only point where the scale factor vanishes, so does S by (|2.18|) .^^ 



Now we consider the condition ii above. In order to maintain a finite Tc, from 
( p.l7|) one should require that the metric function near the singularity should behave 
linearly: 

A{r) -^ -Aoo r + ■ ■ ■ , r ^ oo (2.20) 

where A^o is a constant and the ellipsis may involve (powers of) log's or terms that are 
sub-dominant as r — ;■ oo. This singles out the deep- interior geometry of the type C 



-'^^We ignore the SF term in (2.9) that is subleading in N. 

"'^^This is clear when solving (2.13) by specifying the boundary conditions at r^. See section 3.2 



for more detail. 

"'^^This is the only technical requirement in this section. We note that asymptotically AdS ge- 
ometries are of this type. 

i6\Ye refer to 0] for a derivation of the fact that the BH(TG) metric functions ^(^o) become 
asymptotically the same at the singular points. Thus requirement that exp Aq vanishes is equivalent 
to the requirement that exp A vanishes. 
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above. Then, the rest of the condition i is satisfied automatically: The fact that AF 
vanishes at Tc = oo follows immediately from the definition F{Tc) = E{Tc)—TcS{Tc). 
The enthalpy piece vanish because the entropy vanishes and T^ is finite. The energy 
also vanish by the fact that in the limit T —^ T^ which corresponds to r^, = oo the 
BH becomes of zero size, hence the ADM mass vanishes^^. 

This is of course in accord with the first law of thermodynamics ( |2.19| ) where we 



tacitly made the assumption that AF vanishes at r^. Here, in addition to this we 
learned that, for a continuous transition, one needs Vc = oo and one should require 
( |2.20|) in the deep-interior. This latter condition can be expressed in a coordinate- 



invariant way, as a condition on the asymptotic form of the scalar potential V{^), as 
follows: The Einstein equation (|2.13|) determines the asymptotic form of the scalar 
field: 




^ir)-> .^—A^r + --- , r^oo (2.21) 



where, as in ( |2.2CI| ) the ellipsis may involve logs. Now, the asymptotic form of the 
potential V is determined by the scalar e.o.m., ( p.l5| ), using ( |2.20| ), and inverting 
( p. 21 ) in favor of $ as. 



r($) -> Voo e'^* (1 + K«6($)) , $ ^ oo (2.22) 

where Voo is a constant. 

The term Vgub in ( |2.22| ) is a sub-leading piece: Vgub — ;■ as $ — >■ oo. This piece 
will be important in determining the order of the phase transition in the next section. 
It is interesting to note that, even if log terms were present in ( p.20| ) and ( |2.21| ), they 



conspire to produce exactly the single exponential (with no log $ corrections in the 
exponential) in the leading asymptotics. 

The only dimensionful parameter in the problem is V^o in ( p.22| ). Indeed, V 



has dimensions of £~^ where £ is a typical radius of curvature of the TG geometry, 
for example , it is the AdS radius when the TG geometry is asymptotically AdS. 
Thus, all of the other dimensionful quantities above, i.e. A^o, $oo and T^ should be 
determined in terms of Vr^ 



oo- 



Tc oc Aoo oc $00 oc r^ oc ^/V^. (2.23) 

We shall determine the proportionality constants precisely in the next section. 

To summarize: In the Einstein-scalar set-up, continuous phase transitions of 



order > 2 arise with and only with the asymptotics ( 2.22 ). The transition happens 
at Th = Tc = oo where the entropy, the energy and the free energy vanishes, as 
the black hole has vanishing classical horizon in this limit^®. However, the surface 



''See for a careful computation of the ADM mass. 



17c 

-'^^As explained in section 2.1 above, vanishing of the BH entropy and free energy density is 



equivalent to the vanishing of the entropy and free energy difference between the BH and the TG 
backgrounds, to leading order in l/N"^ . 
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gravity, hence the temperature remains finite. Finally, the value of the transition 
temperature is controlled by constant Voo in (|2.23| ). 



A Comment: We have a BH geometry that has vanishing mass at Tc, to leading 
order in the 1/A^^ expansion. Yet the temperature Tc is finite. One should ask how 
a massless object can have finite temperature? Does it radiate? In fact our situation 
is very similar to what happens in the massless limit of linear dilaton black-holes 



13]. As the authors of [13| show in the special case of a massless linear dilaton 
BH, although the temperature remains finite, the radiation from the event horizon 
vanishes in the massless limit, hence there is no contradiction. In our case as well, 
one can show that the radiation vanishes in the massless limit r/^ — )■ oo. We will show 
this in section § by computing the shear and the bulk viscosity coefficients r] and ( 
and showing that they vanish in the limit r/^ — )■ oo. However, we also show that the 
non-trivial and interesting quantities are not the bare viscosity coefficients rj and ( 
but the normalized ones rj/S and (^/S which remain finite in the limit T — )■ T^. We 
refer the reader to section ^ for more details. 

2.3 Linear dilaton background in NCST 

Now, we consider the special case when the scalar $ is the dilaton of {d + 1) non- 
critical string theory. The NS-NS action of NCST in rf + 1 dimensions read, 

Soc f d'+'xy/^se-'" (i^s + 4(9<l>)2 + ^ + . . .\ (2.24) 

where we only show the relevant terms in the two-derivative sector, the ellipsis refer 
to higher order terms in the derivative expansion. 5c refers to the deficit central 
charge in d + 1 dimensions, e.g. it is 9 — d in a parent lOD critical string theory. 
The NS-NS two- form B is chosen to constant, hence it does not contribute to the 
equations of motion. The subscript s refers to string frame quantities. 

One passes to the Einstein frame by the Weyl transformation g^" = e^^^^'^^^^g^'^. 
The transformation produces an Einstein-frame action of the form ( |2.1| ) with the 
normalization ^ fixed as, 

e = ^. (2^25) 



In this case the requirement for weak phase transitions (|2.22|) become 



4 



V ^V^ e— * + ■ ■ • (2.26) 

Interestingly, the form of the dilaton potential (in the Einstein frame), ( p.22| ) that 
follows from the condition for existence of continuous phase transitions, precisely 
coincides with the dilaton potential that one obtains in the NCST after passing to 
the Einstein frame. At the same time, we know that there is an a'-exact solution to 



the NCST, that is the linear dilaton background. [M. Indeed the asymptotic form 
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of the solutions we found in the previous section, both for the TG ( p.2|) and the BH 
( p.3|) backgrounds^^, become hnear dilaton geometries in the hmit r/^ — )■ oo, if one 
passes to the string frame: 

dsl = [dr' + dx'^ (1 + ■ ■ ■) , ^{r) = ^Aoo r(l + ■ ■ ■ ) (2.27) 

where the corrections denoted by the elhpsis vanish as r — )■ oo. The precise form of 
these corrections depend on Vsub in ( |2.22|) . 

The crucial fact is that when $ is the dilaton of NCST, the string-frame scale 
factor As = A + 2<l>/((i — 1) vanishes when V has the form ( f^.2(j| ), thus one obtains a 
linear-dilaton geometry in the deep-interior, in the string frame. Because the linear- 
dilaton is an exact solution to non-critical string theory in all orders in a' we expect 
that the leading term in the condition ( |2.22 ) for presence of continuous transitions 



is a'-exact. We have more to say on this in section 0. 

In passing, we note that in this particular case of the D-dim. NCST (with 
^ = A/d — 1), the coefficients A^, $oo in ( 2.2CI|) and ( p. 21 ), and the transition 



temperature Tc can be expressed in terms of the parameter V"oo in a precise way (see 
App. for derivation): 

^oo = T^$oo = 4^ % = ^. (2.28) 

~ d-l d-l d-l ^ ' 

This last equation, in particular, grants us with direct control over the transition 
temperature by adjusting the amplitude of the scalar potential in the asymptotic 
region. Clearly, this is very much desired in phenomenological applications. 

2.4 The problem of singularity 

We showed above that existence of a continuous Hawking-Page transition in the 
Einstein-scalar theory requires a singular geometry with a curvature singularity in 
the deep interior. The transition happens at the point when the horizon marginally 
traps the singularity. 

Of course, the Einstein's theory breaks down as one approaches this singular 
limit and one can no longer control the calculations. In fact, in the string frame we 
do not have to worry about the curvature corrections precisely around the transition 
region because the curvature invariants can be shown to vanish in this limit, see 
Appendix |^. Instead the issue is that the dilaton becomes very large hence the 
effective string coupling becomes strong in this region and our assumption of ignoring 
the gravitational interactions cannot be justified. We can avoid this problem, if we 
have the means of fine-tuning the asymptotic value of e* near the boundary to be 
parametrically small. This is what we shall do in the following: we will adjust the 
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In the BH case the blackness function fir) approaches 1 as the horizon size becomes smaher 
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dilaton potential V in such a way that the approach to the singularity is slow enough. 
More precisely, we will choose the subleading terms in the potential such that the 
vicinity around T^ corresponds is away from the large $ region. How to manage this, 



is explained in section O. In section 6 we construct explicit examples where we 
can show that this fine-tuning is equivalent to choosing the asymptotic value of the 
dilaton to be small. 

3. Classification of continuous transitions 

3.1 General results 

In this section we investigate the role of the sub-leading asymptotics in the scalar 
potential Vsubi^) in determining the type of the phase transition. We shall not be 
completely general, but only focus on two interesting cases where Vsub is either an 
exponential or a power-law. 

We state the results before detailing their derivation. We consider two different 
classes of asymptotics: 

(3.1) 
(3.2) 



Case i : Knfe = C e"'^*, k > 0, 


$ -)■ oo 


Case ii : Knfe = C $"", a > 0, 


$ -> oo 


Furthermore, we define the normalized temperature. 




,_T-T, 





T, 



(3.3) 



c 



The asymptotics ( p.l|) and ( p.2| ) the free energy density difference. 

Order n transition : AF{t) ~ t", for k = ^^ ~ \ n>2 (3.4) 

n — 1 

Infinite order transition : AF(t) ~ e~ * " , for a > 0. (3-5) 

In particular when a = 2 one obtains a BKT type scaling in the case (|3.2| ): 

AF(t)~e-^*■^ a = 2. (3.6) 

Another special case is when the scalar field is the dilaton in NCST where the 
normalization of the kinetic term in ( p.l|) becomes C, = 4:/{d — 1). Then, in the 
case ( p. ID above, one has an nth order transition for k = 2/{n — 1). 



3.2 Derivation in a coordinate invariant approach 



We would like to derive the results (|3.4]) and (|3.5| ) in an approach that is independent 
of the chosen coordinate system, where one only needs to specify the asymptotics 
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of the scalar potential as in ( |3.1|) and ( p.2|) . In this section we introduce such an 
approach. The derivation of the thermodynamic properties based on this approach 
is presented in the next section. 

Such a coordinate independent approach was introduced in terms of the so- 
called "scalar variables", in section 7 of 0^*^. The idea is to introduce variables that 
transform as scalars under the diffeomorphism that involve r i.e. r — )■ r{r). Then, 
one reduces the complicated Einstein-scalar system to a system of coupled first order 
equations for these scalar variables. 

One can always show that all physical results, including the thermodynamic 
functions etc. will only depend on this reduced system. 

In case of a neutral black-hole in c? -|- 1 dimensions, one only needs two such 
variables, that we call X and Y: 

xm . 2|, y(*) . i^. (3.7) 

where the constant 7 is given by. 



These functions satisfy the following system of coupled first order equations: 
dX _ ,^ ^2 , .^N /^i , 11 dlogV 



-^-,il-X'^Y)[l^---^), (3.9) 

§^-0(l-A-^.y)I. (3.10) 

We leave the derivation of these and other equations below to App. A. These equa- 
tions are solved by imposing boundary conditions at the horizon: Let us denote the 
value of the scalar at the horizon as ^h- Then regularity of the horizon fixes the 
values of X($/j) and Y{^h) completely in terms of the scalar potential: 

^('^) = $7r$ + ^1 + ^(^'^ - '^)' (3.11) 

X($) = Xh + Xi(<l>, - /) + 0($, - $)2, (3.12) 

with 

(3.13) 
All other sub-leading terms are determined similarly by perturbative expansion near 



^°See also the "thermal superpotential" method in m that is equivalent to the scalar variables. 
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Having determined the functions X and Y by solving these eqs. the metric 
functions are obtained by the definitions (IsTTI) as: 



A($) =Ao + ^ [ ^rf<l, (3.14) 



/($) = exp 



7 / -c?$ 

. J'S'o ^ 



(3.15) 



where $o is the boundary value of the scalar field, and in the second equation we 
used the requirement that the metric function / is unity at the boundary. 

This reduction corresponds to setting $ as the radial variable, with the hidden 
assumption that r($) is a single valued function, i.e. $(r) is monotonically increas- 
ing. 

Finally, one can invert ( p.9| ) to express the scalar potential in terms of the scalar 
variables as, 

Vi^) = '^^^ [l + Y- X^) e-2^/-oo(^W-a^K. (3.16) 

Derivation of thermodynamic relations in terms of the scalar variables is pre- 
sented in App. |B[ Given the scalar potential V{(^) one solves for X($) in (|3.9| ) 
and ( p.lO| ) obtains the thermodynamics of the corresponding BH solution from the 
following two equations. 

AFi^,) = -^ d^^ ^i^-mi.'^.) ^, (3.17) 

4<-^D y$„ d^h 

n^u) = , J ,, e^(*'^) V{^,) e^ /*>(*) ^*. (3.18) 

47r(a — 1) 



The first one follows from the 1st law of thermodynamics ( |2.y| ) and the expression 
for the entropy ( p.llj ) . The second one follows from expressing ( p.4|) in terms of the 
scalar variables. We refer to Appendix ^ for details. 

The constant I that appears in these equations is the length scale, (^.23]), that is 



given by the curvature radius of the geometry at the boundary $0) which is in turn 
determined by the value of the potential at $0^ 

For example, in case of pure AdS, the scalar should be constant, hence X = and 
from (p.l9|) follows the standard expression for the cosmological constant in d+1 
dimensions. 
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4. Scaling of thermodynamic functions 

4.1 Free energy 

Now, it is a matter of straightforward calculation to see how the free energy scales 



with temperature for asymptotics of the type ( ^.11 ) and ( |3.2D that hold in the region 
1 ^ $ < $/i. One first solves X and Y from (|3.9|) and ( |3.10 ) and shows that, 



X{<^) = --^ + 6X{<^), 



Y{'^) = Yo{<^) + 6Y{<^), 



(4.1) 



where SX and 6Y are sub-leading with respect to the first terms in the limit 1 ^ 
$ < $fe. 

It turns out that we do not need the precise expressions for 6X, Yq and 6Y above 
to determine the type of the phase transition. However they become important, for 
example, in determining the viscosity coefficients around T^. We present them, in 
Appendix 0. 

One determines the scale factor A($) from ( p.l4| ) as, 

(5X($) d^ + const. (4.2) 



7 



The constant piece is unimportant in what follows and one can show that the second 
term is sub-leading with respect to the first one, see Appendix ^ 

Substituting ( |4.1| ) and ( [4.2[ ) in ( p. 181 ), one finds that all the dependence on SX 
cancels out and the rate of approach of T toward T^ is directly determined by Vgub- 



Using the definition of the normalized temperature ( p.3|) , one finds: 

t = Vsubi^h). (4.3) 

Now the free-energy as a function of t is given by ( p.l7|) and it can be rewritten as. 



AF(t)oc / rfte(^-i)^®. (4.4) 

Jo 

The dependence of the scale factor on t should be found by inverting ( [4. 3] ) and 

substituting in (the leading term of) (|4.2|): In the cases ( |3.1| ) and ( |3.2| ) one finds, 

1 



Case i : 



Ait) 



K 



e 



d- 1 



Case ii : A{t) 
The free energy follows from ( 



e 



d-1 



hgit/C) + ■ 



(t/cy 



t ^0+ 



t^o^ 



(4.5) 
(4.6) 



as 



Case i : AF(t) oc t'^^^+\ 
Case ii : AF(t) oc 6^'*"*^^+^, 






(4.7) 
(4.8) 
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where C = \/^{d — 1)C^ in the second equation. We see that AF vanishes, as it 
should, for arbitrary but positive constants ^, k and a. 
In the special case of 

K = y/ad-l)/{n - 1) (4.9) 

in ( |4.7] ), one obtains an nth order transition in case i. Moreover, if the scalar field is 
dilaton, then ^ is given by ( p. 251) and (|4.9|) reduces to. 



n — 1 
In this case the required asymptotics in ( p.22| ) for an order-n transition becomes 



(4.10) 



V^($) -> e^^* + Ce("-i)(''-i), $ ^ cx). (4.11) 

For example, for a 3rd order phase transition in a dual 2+1 dimensional theory, one 
should require, 

V{<l>) -> e^* + Ce", $ ^ oo. (4.12) 



On the other hand, the special case of a = 2 in ( p.2| ) is a BKT type transition 
where the free energy scales like^^ 

-c' 
F{t) ~ e ^ . 

We note that the setting is rich enough to cover different BKT type scalings that one 
encounters in the condensed- matter literature for different values of a. For example. 



in the dislocation theory of melting in two-dimensional solids [^ one encounters a 
different value, a = 2.70541 ... in (|4^ ). 

We also note that this set-up covers other, more general type of divergent tran- 
sitions where one of the derivatives of the free energy becomes not discontinuous, 
but diverges. For a recent reference, cf. |jl6[. Clearly, these type of transitions are 
ubiquitous in the case i above, for arbitrary real k. 

4.2 Fine-tuning the potential 

Before continuing further, we would like to comment on the issue of the transition 



region coinciding with the singular region in Einstein's gravity, section |2.4| . In this 
paper we are only interested in how the thermodynamic quantities scale with the 
reduced temperature t and not in the overall size of these quantities. In this case, 
one way to avoid the problem of running into the singular region is to fine-tune 
the parameters of the potential V{^) to parametrically separate the singular region 
<I>/i ^ 1 from the vicinity of temperatures where the scaling sets in, i.e. t <C 1. 
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One can also get rid of the power-like corrections in (4.8) by slightly modifying the asymptotics 



in (p.2[) by adding log corrections, but we shall not consider this here. 
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Equation ( |4.3|) suggests an obvious way to manage this separation just by choos- 
ing the phenonienological parameter C in ( p.l|) and ( p.2|) , small enough. In the 
Einstein frame how small it should be chosen is determined by the demand that the 
size of the curvature R(? is not too large^^. Suppose we decide to trust Einstein's 
gravity up to the point where the curvature becomes some order 1 number A. Then, 
in the region where V is approximately given by ( [^.22| ), one has, 

R(xV^ e^^*'' ~ A/i"^, (4.13) 

where the proportionality constants are 0{1). We refer the reader to Appendix ^ for 
a derivation of the curvature invariants. Substituting this in ( |3.1| ) or ( |3.2| ) determines 
how small the constant C should be chosen. In case 1, for example, one needs. 

In our case C is a phenomenological parameter that can be chosen arbitrarily small. 

Another means to accomplish the separation of the transition region t k. {) and 
the singular region $/j 3> 1 is to tune the boundary value of the dilaton to be very 
small. This can be done in an asymptotically AdS solution, see section 6. 

Either of these ways seem to be a crude way to avoid the problem and there is 
a prize that one pays. For example when we consider quantities given by ratios such 
as 7]/ S and C,/ S, we will see that there is an ambiguity in determining their value at 
T 

4.3 Entropy and specific heat 

As we focus on continuous transitions in this paper, the classical entropy difference 
vanishes exactly at the transition^^. The scaling of the entropy near Tc is determined 
by dmi) and (pj), | 



Casei: /\S{t) (x t"^^^ + 0{N-^), t ^ 0+ (4.15) 



Case ii : AS{t) oc e"^ * " + O^N''), t -^ 0+, (4.16) 



where C is defined below ( [4.8|) . 



^^One can address the same issue in the string frame. In this frame the curvature invariants vanish 
in the transition hmit (see Appendix M, however a large vahie of the dilaton signals breakdown 
of our assumption of weak gravitational interactions. In section 6, we construct explicit examples 
and show that the condition derived here is essentially equivalent to choosing the boundary value 
of the dilaton small enough so that it stays small in the transition region. 



^''As discussed in section 2.1, by vanishing of the entropy we really mean it becomes of 0{N ^) 
in the large N limit. 
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The specific heat is obtained from ACh = TdAS/dT as, 

Casei: AC^{t) oc t"^^^'^ + 0{N-^), t ^ 0+ (4.17) 

Case ii : AC^{t) oc e-^'''^r^'^ + 0{N-^), t -^ 0+, (4.18) 

Clearly, for any k > \/^{d — 1) the specific heat diverges at the transition. In the 
special case of (|4.10|) with n = 2, it becomes discontinuous, as it should in a second 
order transition. 

4.4 Speed of sound 

The (normal) speed of sound is defined as c^ = dP/dE. In isotropic systems that 
we consider here, the pressure is given by the free energy density and the speed of 
sound in the BH phase is most easily determined from 

cIbh = ^- (4.19) 

Now we ask the question, what does c^ ^^j become in the limit T — )■ Tc"! Equation 
( [4.19| ) can be written as. 



2 AS + Stg 



For the case of a 3rd or higher order transition, the leading order terms AS and 
ACj, vanish as in section [4.3| and one is left with the terms of 0{N~'^) both in the 



numerator and the denominator. Their ratio is 0{1) and it is given by Stg/Cv,tg 
which is nothing but the speed of sound in the thermal gas phase. Therefore the 
speed of sound stays continuous across the transition, which is indeed what should 
happen for a transition of order higher than 2. In the case of a second order transition 
A^* vanishes, but AC^, stays constant and one finds a discontinuity in the speed of 
sound. All of this is completely in accord with a system that obeys the laws of 
thermodynamics . However there is an important ambiguity which we discuss below. 

4.5 An order of limits issue and divergence in the dilaton 

An ambiguity arises when we consider the A^ — ?■ oo limit in the discussion above: We 
observed that taking the t — )■ limit while keeping A^ finite one obtains cj. ~ 0{1) at 
the transition point. Now, one can further take the N ^ oo limit and c^ still stays 
the same. On the other hand, had we taken the A^ — )■ oo limit first, then the 0{N~'^) 
terms in ( [4.20| ) would drop out and we would obtain the result c^ ~ t around the 
transition. Thus, we learn that the limits A^ — )■ oo and t — )■ do not commute. 

The resolution of the puzzle is clear: one cannot really take the t — )■ limit for 
finite A^ without taking into account the gravitational corrections. This is because 
the Einstein's theory breaks down in this limit. For definiteness, let us consider 
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Einstein's gravity as an effective theory to a d + 1 dimensional non-critical string 
theory^^. Then, the effective Newton's constant at point r is given by 

Geff = (167rMp^-^Ar2)-ie2*('-) ~ (?2^f-^e2*W, (4.21) 

where Qs is the value of the dilaton field on the boundary. To be definite, we con- 
sider an asymptotically AdS geometry where the dilaton approaches a constant value 
$(r) — )■ $0 as r — )■ and Qs = exp($o)- This implies that, for finite N there is a limit 
in r/i beyond which we cannot ignore the gravitational interactions. In the specific 
case of ( p.l|) it is determined from ([4.21|) "^^^ 



as 



tum ~ C(4v^(M/)('^-^)/2^)-^ (4.22) 



where C is the constant defined in ( |3.1| ). Beyond this limit one should take into 
account the gravitational corrections. 

One can ask whether one can send tum to zero by fine-tuning a parameter in 
the theory. One can indeed fine-tune tum to be arbitrarily small by choosing C or 
1/N^ ~ e*° to be small. Sending C to zero is not desirable, as in this case the 



subleading term in (|2.22|) vanishes and the resulting potential displays no transition. 
But one can send e*" to zero instead. This corresponds to the strict N ^ oo limit. 
The conclusion is that for a system that corresponds to finite N, regardless how 
large N is, one should consider 0{Gd) corrections in order to calculate the value 
of a quantity such as c^. On the other hand, the expected result that c^ becomes 
discontinuous at T^ in a second order and continuous in a higher order transition 
should still follow after taking into account the 0{G£,) corrections because AS and 
AGc should still vanish at Tc also after the corrections. The only ambiguity is in 
determining the value of a quantity given by a ratio such as c^, precisely at Tc. We 
encounter other examples of such quantities below. 

5. Viscosities and dispersion 

Viscosity coefficients are important observables probing the dissipative behavior in 
the dual field theory at finite temperature. In our case there are only two viscosity 
coefficients, that determine the response of the system to shear and bulk deformations 
respectively. 

5.1 Shear deformations 

The shear viscosity is given on the gravity side by the flux of the "shear" gravitons, 
i.e. the traceless-transverse component Htt of the metric fluctuations. Due to the 



^''The precise nature of this theory is irrelevant for the sake of discussion. One could also consider 
a compactification of II(A)B on a 9 — d dimensional internal manifold. 
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we considered the dimensionless effective Newton's constant Ge//^'' ^■ 
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well-established universality in the case of two-derivtive gravity [17|[1S] the shear 
viscosity to entropy ratio in the BH phase is given by, 

where the ellipsis refer to a' and 1/N corrections. We note that at the classical 
level in the two-derivative Einstein-scalar gravity no temperature dependence arises 
in this ratio, and the result is simply given by (|5.1|) for at any T ^ T^. 

What happens exactly at T = Tc is similar to the story of the speed of sound 
above. To be more precise by including the l/N"^ corrections, one should write 
V = Vo ~^ ^V where rjo is the leading order piece in the l/N"^ expansion, 770 = AS/An 
and 67] refers to l/N"^ corrections. One has. 

In the limit T -^ T^ the leading order terms vanish, but one still obtains a 0{N^) 
answer, 

"^ -^^- (5.3) 



Sbh Stg 

In order to determine the value of this, one has to include the (9(1/A^^) corrections in 
the analysis. The situation is different in the strict A^ — )■ 00 limit where the (9(A^~^) 
pieces in (|5.2|) drop out and one retrieves the universal value tj/s = l/An also at the 
transition point. 

5.2 Bulk deformations 

The hulk viscosity on the other hand is more interesting as it displays T dependence 



already at the two- derivative level |jT9|[^. It is given on the gravity side, by the 



flux of the rotationally-invariant fluctuations of the metric h^ = hu = /i22 = ■ ■ ■ = 
hct-i^d-i- The computation is much more involved than the shear case, due to mixing 
of these fluctuations with the fluctuations of the scalar fleld 5(j). An ingenious method 



to disentangle this mixing is devised by Gubser et al in [0. The main idea is to 

consider a gauge where the fluctuations of the scalar is set to zero. This is equivalent 

to carrying out the computation by using as the radial variable. The method is 

described in detail in [|l^ for the case of 0? = 4. Here, we generalize this calculation 

to arbitrary dimensions. This is straightforward, yet technically involved and we 

present the details of this computation in App. ^ . We shall only highlight the main 

steps and results here. 

Kubo's linear response theory |^ relate the bulk viscosity to the retarded Green 

function of the stress tensor as ^^ 

4 1 

^0 = ~ 7a — T\2 ^^"^n ~^"^ Gr{u) (5.4) 
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The subscript refers to the leading term of the quantity in the 1/iV^ expansion. 
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where the Green function is evaluated at zero momentum and it is defined by, 



Gr{uj) = -i / dtd''-'xe'^'e{t){[-Tl{t,x),-THt,x)]). 



(5.5) 



and i,j denote the spatial components only. The imaginary part of the Green func- 
tion maps to the fiux of the hb gravitons on the gravity side. It is determined by a 
long and non-trivial computation (cf. [0 and App. ^ of this paper). The compu- 
tation is carried out in the gauge 5$ = 0. In other words $ is taken as the radial 
variable. We denote the metric functions in this gauge as, 

ds^ = e2^(*) (-/($)rft2 + dxl_,) + e2^(*)-'^'^' 



fm' 



The final result reads. 
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where prime denotes d/d^. The metric fiuctuation h^ satisfies the following second 
order equation : 

,2 // £ f 

f 



^'^^'' {-^^ + ^^ " ^)^' " ^^' + 7) +^ 



, 1 ^e'^-'^ + Cb' 



d-lfA' 



0. 



(5.8) 

Here prime denotes d/d^. In passing we note that this equation looks simpler and 
more transparent in the r- variable: 



/., + /., I (ci-l)i+^ + 2|j+/^J^-^| 



0. 



(5.9) 



where dot denotes d/dr now. We note that reference to the normalization of the 
kinetic term ^ disappears in this form. We also note that it reduces to the shear 
fiuctuation equation for X = const. 

Equation (^.9|) is solved by the boundary conditions hj, = 1 on the boundary and 
the in-falling boundary condition at the horizon: 



hb -> Cb(rh-r) *4^T. 



r -^ rh- 



Here Cb is a T-dependent multiplicative factor. Inserting ( ^.101 ) in 
variables to $) yield, 

ImGRiuj) = 






Now, the bulk viscosity follows from the Kubo's formula ( p.4| 

Sbh Svre'''" \V^) 



(5.10) 
(and changing 

(5.11) 
(5.12) 
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5.3 Universality at Tc 

In the limit $/i — t- oo i.e. T ^ Tc one obtains a universal result for the bulk-viscosity- 



to- entropy ratio. As shown in pO[ the coefficient |c(,| approaches to 1 in that limit. 
The ratio V'/V also becomes a constant that is given by the exponent in equation 
( g:22|) . Thus, (|5A2|) becomes, 



Co 



S. 



BH 



-, (5.13) 



27r(rf-l)' 

at the transition. We note that this is a universal result in the sense that it depends 
neither on the value of Tc nor on the details of the underlying field theory system, 
nor on the order of the phase transition. The details of the computation that leads 
to this result are presented in App |^. 

As for the shear case, the result is valid only in the limit where one can ignore 
gravitational interactions. Just as in the discussion below equation ( |5.1| ) above, if 
there exists 1/A^^ corrections then the result changes not by 0{N~^) but by 0{N^). 
Denoting the l/N"^ correction to Co as ^C, one has, 

C ^ Co + '^C _^ _^ as T ^T (5 U) 

Sbh as + Stg Stg 

The latter is determined by computing the l/N"^ corrections in the background and 
there is no reason that it should coincide with ( [5.13| ), even in the large N limit. 
However if one considers a system with N = oo strictly, then one should take the 
N -^ oo before t — )■ and the result can be trusted in principle. In the next section 
we construct explicit kink solutions that flow from an AdS on the boundary (with 
a fixed value of e*°) and linear-dilaton in the deep interior. In these solutions one 
can indeed send e*° — )■ in a consistent fashion and one obtains the universal result 
above. It would be extremely interesting to find a similar solution in IIB and check 
whether the result given by ( |5.13| ) holds there. 

5.4 Scaling of viscosities near Tc 

There is no ambiguity however, in the scaling of the viscosities rj and ( near Tc, as 
this is determined by the leading piece in the 1/A^^ expansion. From the results 
above, we find that both rj and ( scales as the entropy density near Tc. 

7], Coc A^(t)oce('='-^)^W. (5.15) 

The latter is given by (|4.5|) and (WM) in cases i and ii: 



Case i: r], ( oc t « t-J-O (5.16) 

Caseii: r^, C oc e^'*"", t ^ 0, (5.17) 

where C is a constant given below equation ( [4.81) . The latter is the typical scaling 
of the order parameters in the Kosterlitz-Thouless transition |ID], [P^ . 
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5.5 Dispersion relation of pressure waves 

Now, we are in a position to write down the dispersion relation of pressure waves 
near criticality, including dissipative terms. The hydrodynamics predict a dispersion 
relation of the form, 

i f d r/ , 1 C \ 2 



T \d- ISbh 2SbhJ 

From the results in the previous two subsections, eqs. (|5.1|) and (|5.13|) , in the strict 
N = oo limit one finds, 

Casei: ^ ^ ^' ^~ ^^ [d^) ^^' ^ "^ ° ^^'^^^ 

^ i f d + 1\ r, , ^ 

Caseii: ^ ^ ^" ^ ^ 4^^ ( 73^ ) ^ ' i ^ 0, (5-20) 

where we only show the leading terms in t in both of the terms. One should remember 
the ambiguities in these formula for finite N as explained in the previous subsection. 

6. UV geometry and analytic kink solutions 

So far we have not specified the geometry near the boundary. The only assumption 
we made in the derivation of section 2 is that A'{r) < near r ^ 0. This includes the 
case of asymptotically AdS backgrounds. In fact, for holographic applications this is 
the most desirable case. Therefore it would be interesting to find explicit solutions 
that fiow from an AdS geometry in the UV to a geometry of the type given by the 
solution of the IR potential ( |2.22| ) with Vsub given by either of (|3.1|) or ( |3.2| ). In this 
section we construct such explicit examples for the case of ( |3.1| ). This will the most 
interesting case for applications to holographic superfiuids in 0. 

An asymptotic AdS geometry is also desirable for another reason. In the case 
when the scalar can be regarded as the dilaton of NCST, its value goes over to 
a constant $(0) = $0 in an asymptotic AdS geometry which provides a tunable 
parameter. On the other hand, the effective Newton's constant is related to this 
value as, 

Gd oc e2*o. (6.1) 

Therefore, this value can be tuned to make the gravitational interactions very small. 
In other words, $0 wiU be the parameter that corresponds to the "color" N in a dual 
field theory 

N^ oc e-2*o. (6.2) 

The large N limit corresponds to tuning $0 to be very large. 

At the same time, we argued in section ^^ that in order to keep the transition 
region away from the curvature singularity, one needs another fine tuning C ^ 1 
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where C is the constant defined in (3J-). In principle, these two separate fine-tunings 
may contradict each other. Here, we will show by constructing explicit examples that, 
on the contrary, these two fine-tunings are essentially the same^^. We shall focus on 
the case where the scalar field is regarded as the dilaton and its normalization in 
(Oisgivenbye = 4/(d-l). 



6.1 Two-exponent potentials 

For notational simplicity we define the function A 

A = e*, Ao = e*o. (6.3) 

Consider the potential, 

V{\) = V^\^^ {l + C\-^). (6.4) 

According to the discussion in section |^, the black-hole solution that follows from this 
potential exhibits a continuous type Hawking-Page transition with the order of the 
transition determined by the constant k as in ( |4.1CI| ). For example the transition is 



second order for k = 2. The transition region is in the asymptotic IR region A/^ S> 1. 
We observe that this simple potential allows for an AdS extremum at^^ 

A„ = ci(^i<^-ll" (6.5) 



when K[d — 1) > 4. As a result, the two-exponent potential ( \6.4\ ) should have a kink 
solution that flows from an asymptotically AdS geometry in the UV at X = Xq towards 
the asymptotically linear- dilaton geometry in the IR at X = oo. 



We observe from (|6.5|) that the large N limit (|6.2|) corresponds to C ^ 1. This 
is reassuring because the same limit is required in section ^^ in order to separate 
the large curvature region from the transition region. 

One can construct this solution analytically for the thermal gas case ( p.2|) , for 
special values of k and d. This can easily be done by the method of scalar variables. 
Appendix ^^ as follows: The TG solution corresponds to the special solution of 
( P.10|) Y = 0. In this case the Einstein-dilaton system reduces to a single first order 
differential equation, 

dXo _ ,, ^2N /^. , 1 1 dlogV 



« ^<^-'^«n'+^^^^'' <'■'' 



^^This is not hard to see directly from (3.1) where shifting the vev of $ can be regarded as tuning 
the constant C. 

^*We assume C > 0. 

^^It can also be obtained by the method of "fake super-potential", [ESl. 
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where 7 is given by ( |3.8|) , 

A very nice coincidence is that one can construct an analytic solution to (|6.6|) pre- 
cisely in the most interesting case, namely k, = 2, d—1 = 3. This case corresponds to 
a second order phase transitions in three spatial dimensions, indeed the most relevant 
case for many condensed matter applications! 0. The solution reads. 

In the UV limit Xq — )■ which indeed corresponds to an asymptotically AdS ge- 
ometry^°. To see this, one recalls the definition of Xq in ( p.7|) and vanishing X 
corresponds to constant dilaton. 

One can construct the metric and the dilaton in the r-frame explicitly by solving 
equations (|A.5|) and ( |A.6|) with X = Xq given by (|6.8|) and Y = 0. The result is. 



4 ch3f^) 
sh (§) 

A(r) = Ao ch(-). (6.10) 

Clearly, the geometry becomes AdS near r = 0. The metric in the string-frame is 
obtained from ( |6.9| ) by the Weyl transformation ds^^ = X^'^ds^: 

dsl = coth2(^) K + dxl, + dr') . (6.11) 

In the IR region r — t- 00 the background in the string-frame becomes linear-dilaton: 

dsl^de + dxl_, + dr\ *W^^- (6.12) 



Therefore the background ( p.ll| ) and ( |6.1CI| ) is a solution to the 5D Einstein-dilaton 



theory that interpolates between an AdS background near the boundary and a linear- 
dilaton background in the deep interior in the string frame. 

It would be very useful to obtain an analytic black-hole solution to the potential 
( |6.4| ) by turning on a blackness function in ( [6.91) . In the language of scalar vari- 
ables ( |3.9| ), ( ^.101 ), this corresponds to turning on a non-trivial Y . Now, the system 
becomes more complicated and unfortunately we could not find an analytic solu- 
tion. However a numerical solution can easily be constructed by solving the system 
by a numerical code. We checked that the black-hole obtained this way obeys the 



general conditions described in section 2.2, in particular it asymptotes the TG kink 



above in the limit r/^ — )• 00. To our knowledge, the closest analytic BH solutions are 

the Chamblin-Reall solutions to single exponential potentials [^. However, these 

solutions do not exhibit any Hawking-Page transitions. 

■^"Equation ( |6.6| ) looks singular in the limit Xq — > 0. In fact it is not as the pole of Xq at Aq is 
cancelled against the pole in V'{\) 
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6.2 Three-exponent potentials 



The analytic solution found above is for the particular case n = 2, d — 1 = 3. It is 
interesting to search for similar analytic solutions for an arbitrary x and d. This is 
not possible with a simple two-exponential dilaton potential. Instead we consider a 
three-exponential potential of the form: 



vi\) = VocX~ 1 + ex-"" + ex 



-2k 



(6.13) 



This potential possesses AdS extrema for generic values of k and d. Let us again 
consider a thermal-gas type ^]2| ansatz. One can find an analytic solution for the 
special value of C: 

In fact this class of potentials can be regarded as a generalization of ( p.4p with k = 2 
and d — 1 = 3. In this special case the constant C in (|6.14| ) vanishes and the three- 
exponential potential ( |6.13D reduces to i 
reads, 



vd\-+w 



Solution of (|6.6|) in the general case 
1)-1. (6.15) 



1^" 



One checks that, this solution reduces to ( |6.8| ) in the special case oi ht = 2, d — 1 = 3. 
The special point Xq corresponds to the extremum of (|6.13| ) and it is given by. 



An 



Cq 



K 



(6.16) 



The solution ( |6.15| ) corresponds to a kink that interpolates between an AdS geometry 
at Ao to an asymptotically linear-dilaton background at A — )■ oo. 

We again observe that by tuning the constant C to small values one can establish 
two goals at the same time, namely to separate the transition region t fa and the 
high curvature region Xh ^ 1 and to make the effective gravitational coupling weak. 

It is not possible to obtain a closed form expression for the metric in the r-frame 
31 for generic k, and d but it is easy to write down the metric directly in the A-frame: 



ds^ = X-T^il-^^'' 



dt^ + dx'^^_^ + Cq 



A'^ 



A-2 dX^ 



(6.17) 



where q is given in terms of k and d in ( |6.15| ) and the constant Cq is given by, 

2£ ^ ^ 



q{d-l) 



(6.18) 



31 



At best one can express r in terms of A as a combination of hypergeometric functions. 
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The metric is defined in the range A G (Ao,C)o). It is easy to see that the metric 
becomes AdS near A ~ Aq by a change of variable A — Aq oc r'^ near r ^ 0. On the 
other hand, for large values of A it asymptotes to the linear-dilaton background in 
the string frame: 

dsl ^ dt' + dxU + (^-^^^^^ rf$l (6.19) 

The solution given by ( |6.15| ) constitutes a one-parameter family of kink solutions — 
parametrized by k, — that interpolates between AdS and linear-dilaton backgrounds 
in the string-frame. The geometry is regular everywhere in the string frame. It 
would be very interesting to embed this geometry in IIB. One may be able to do so, 
by finding an analogous solution in a consistent truncation of IIB with non-trivial 
dilaton. Apart from our main interest in this work, namely that black-holes on 
these backgrounds exhibit continuous type Hawking-Page transitions with the order 
parametrized by k, they may have other applications in different contexts. 

7. Stringy corrections 

In this paper, we consider two- derivative Einstein-scalar gravity in the limit of small 
gravitational corrections, Gd — > 0. On the other hand, we showed that the continu- 
ous phase change occurs in a singular limit, where the black-hole horizon marginally 
traps a curvature singularity. A crucial issue then, is the possible a' or Gd corrections 
at this point. We comment on these issues in this section. 

First of all, we should stress that the results that are obtained above are con- 
sistent and robust in the framework of Einstein's gravity. The two-derivative theory 
breaks down only near curvature singularities, or when the dilaton field becomes 
very large. On the other hand the curvature singularity is always cloaked by the BH 
horizon as one approaches the transition point t — )■ O"*", and the dilaton can always 
be kept from growing too large as long as t 7^ strictly by fine-tuning the dilaton 
potential. How to manage this is explained in section [4.2| . Therefore the questions, 
as well as the answers that are studied so far are definite within this framework. 

On the other hand, a basic motivation for us to study continuous phase transi- 
tions in gravity concerns possible applications in the dual gauge theories or condensed 
matter systems through the holographic correspondence, and a robust framework to 
relate the two sides of the correspondence is string theory rather than Einstein's 
gravity. Therefore, one should consider possible stringy corrections if one desires to 
interpret the results of this paper from an holographic point of view. 

7.1 Higher derivative corrections 

The stringy higher derivative corrections to Einstein's gravity are measured by the 
typical curvature of the geometry Ri'j. in units of string length ig = y/a'. Concerning 
these corrections, one can ask two separate questions: 
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1. Are the conditions derived in sections p.2| and El for the existence and the type 



of the continuous transitions still valid once we include the cJ corrections? 

2. For which quantities and in which regions can we trust the observables if we 
ignore the a' corrections? 

Let us consider the first question above. For existence of a continuous phase 
transition, clearly, the (a' corrected) geometries should asymptotically become the 
same in the IR region. For a BH geometry that asymptotes a TG geometry near 
the boundary, this condition can most generally expressed as vanishing of the ADM 
mass of the BH. 

Let us now specify to the case when the scalar is the dilaton field of the d + 1 
dimensional non-critical string theory. Then the normalization of the kinetic term is 
given by ( |2.25| ). The required asymptotics for a continuous phase transition is given 
by (12:261) : 

In the string frame Qs^^u = exp(2$/(i — l)g^u the action reads 

^^ = JTTT^ [ d'^-'xV^.e-'" {R, - Ud'^f + K($) + ■■■), (7.1) 

IG-kGd J 

where the subscript s refer to string frame quantities and the ellipsis denote the 
(unknown) higher derivative corrections. Here the string frame potential is related 
to the Einstein frame potential by, 

F($) = K,($)e^*. (7.2) 

Therefore the requirement for existence of a continuous transition now becomes Vg — )■ 



const, in (|7. 1| ) , as $ — t- oo. The remarkable fact is that, the unique solution to this 
asymptotic potential with the TG ansatz (or with the BH ansatz in the limit $/^ ^ 1) 
is nothing else but the linear-dilaton solution, that is a' exact! Therefore, we expect 
that the condition ( p.22| ) be a' exact in this case. 

There is still a possibility of a' corrections in the sub-leading terms ( p.l| ) and 
( p.2| ). This means that the constants k and a in these equations can get renormalized 
by the a' corrections^^. However, clearly this would only affect the definition of 
these constants. In other words, an nth order transition would be given by /« = 
2/(n - 1) + 0{a') instead of (|4lO| ). 



Finally, in the case when the scalar is not the dilaton of a non-critical string. 



then there is a possibility that also the leading term (|2.22| ) be corrected. Avoiding 
this is another benefit of identifying the scalar with the dilaton oi d + 1 dimensional 
non-critical string theory. 



■^^ There is also the possibihty of power-law corrections in (|3.lD and log-corrections in (pj 
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Now we move on to the second question above. First consider the thermodynamic 
functions: The entropy of the black-hole should receive the usual ol corrections. 
These are of two different origin: a) The definition of the entropy as the horizon size 
will be corrected by Wald's generalization |^; b) the background geometry itself will 
receive ol corrections. The definition of the surface gravity, hence the temperature, 
is subject to similar corrections of two separate origin. However, the first law of 
thermodynamics should still be valid, therefore we expect that the equation ( ^.17| ) 
(or (|2.19|) ) should survive the corrections. Clearly, all of the other thermodynamic 
quantities that we obtained in section ^ are subject to a' corrections in the region 
of high curvature in the Einstein frame. We expect these corrections to renormalize 
constants k and a in (|3.1|) and (|3.2| ). 

Secondly, one can ask if the thermodynamic quantities are subject to higher 
derivative corrections away from the transition region. The size of the corrections 
away from the singularity should be determined by the ratio Isl^ where i is the 
typical size of the asjTiiptotic geometry, e.g. it is the AdS radius, if the near boundary 
geometry is asymptotically AdS. This ratio is holographically related to an observable 
in the dual field theory. For instance, it is proportional to the gauge string tension 
in the holographic QCD models 0, whereas it is related to the interaction strength 
between the spins in spin-models 0. Then, the issue of corrections to thermodynamic 
quantities become a phenomenological question, and the answer depends on the 
particular dual field theory that one desire to describe holographically. For example, 
in the holographic QCD models considered in 0, it was found by matching the string 
tension on the lattice-QCD that the ratio is (9(10^^) and the a' corrections are small. 

Finally, we note that the quantities that are naturally evaluated in the string 
frame should be trustable in the vicinity of the transition because the curvature 
invariants in the string frame in the transition region becomes very small as shown 
in Appendix ^. Such quantities involve F-string and D — brane configurations in 
the probe approximation. In the companion paper, 0, we map these configurations 
to spin-correlation functions in a dual spin-model. Therefore, we expect that these 
important observables should be immune to a' corrections — except that corrections 
may arise through dependence on the scaling constants k and a in { \i.l\ j and (|3.2|). 
We refer the reader to for a more detailed discussion. 

7.2 1/A^^ corrections 

Similarly the gravitational interactions for finite Gd also becomes important near 
the transition region t ~ 0. This is because the dilaton field^^ diverges precisely at 
t = 0. This means that the effective string coupling Qs becomes large in this region, 
unless N can be set to oo in a strict manner. 



■^■^One should consider the case when the scalar is not dilaton separately, see the discussion in 
section ||. 
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If in the dual field theory one has finite N — regardless how large it is — our 
calculations can only be trusted in the region 



For a derivation of this bound, we refer to section |4.5| . Beyond this point, one should 
take into account corrections. Therefore, in the case of large but finite A^ (small but 
finite Gd), what one has to do in order to compute the rational quantities at Tc, is 
to compute the spectrum of fluctuations in the TG geometry, and from this compute 
the determinant correction in the saddle-point approximation, in the limit T — )• T^. 

We avoided this problem here by tuning the boundary value of the dilaton e*° oc 
A^~^ to be parametrically small so that the restriction above is pushed towards t — )■ 0. 

On the other hand, we may be interested in the strict N = oo case for a physical 
reason. For example, in we propose that the strict N = oo case corresponds to 
the XY-spin models that constitute canonical models of superfluidity. Then, if we 
are interested in a holographic model for superfluidity, we can take the strict A^ — )■ oo 
limit. In this case one can extend the safe region up to t — )■ O"*". However, there is 
an analogous problem where the UV region becomes highly curved. 

8. Discussion 

There are several possible directions one can improve on the results obtained here: 

• Holographic implications in dual field theories. 

Hawking-Page transitions in gravity generically correspond to confinement- 
deconfinement transitions in the dual gauge theory |^. Therefore, the back- 
grounds presented here are natural candidates for continuous type confinement- 
deconfinement transitions. More interestingly, it is also possible to relate grav- 
ity to field theories that describe continuum limit of certain condensed matter 
systems such as the spin models. We refer the reader to our companion paper 
for a detailed exploration of this latter application. 

• Embedding in string theory 

One of main results is that the continuous transition corresponds to the limit 
where the BH horizon marginally traps a curvature singularity where the Ein- 
stein's gravity becomes not trustable. We avoided this problem by fine-tuning 
the coefficient of the subleading term in the potential to be very small. This 
crude way serves our purpose to determine the qualitative scaling behavior of 
observables in the model, but clearly it is desirable to obtain better control over 
the model for holographic applications. In order to resolve the singularity one 
may think of embedding the background in the full string theory. Indeed we 
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showed that the deep-interior asymptotics of the BH geometry near the transi- 
tion is a hnear-dilaton geometry in the case when the scalar is a dilaton, which 
is an exact solution in the string frame, in all orders in a'. This observation 
should be relevant in embedding the set-up in the full string theory. However 
we need more than this because the full flow geometry is not a' exact and there 
is no obvious way to study this solution in a systematic expansion in the full 
string theory, because the typical curvature scale in the geometry is i/£s ~ 1. 
Therefore, a better strategy is to start with a consistent truncation of a critical 
string theory and directly search for continuous critical phenomena there. A 
helpful starting point for this investigation may be the work of ||26|p^ 

We also observe that the asymptotic form of the scalar potential in ( p.22| ) 
and ( ^.11 ) is sum of exponentials that quite generically appears in consistent 
truncations of IIB and IIA critical string theory. We leave this investigation 
for future work. 

The case when the scalar is not dilaton 

Considering the scalar field as dilaton has various advantages. First of all, if 
one would like to think of the framework as an approximation to non-critical 
string theory^^, then the most economical choice is to consider $ as the dilaton 
field of NCST. Secondly, for this choice the IR asymptotics become the linear 
dilaton geometry that is known as a full solution to NCST. As we show in P, 
another reason is that, in this case scaling of certain observables in the dual 
field theory, such as the magnetization if the dual theory is a spin-model, has 
expected scaling properties as a function of t. Yet, it would be interesting to 
investigate the cases when $ is just a scalar field, and work out the implications 
in a hypothetical dual theory. It is also interesting to search for continuous 
Hawking-Page transitions when there are more than one scalars present. 

The specification of the UV geometry 

In this work, we were mostly interested in scaling of thermodynamic observables 
near the transition that corresponds to the IR limit of the geometry. For the 
various applications in holography however, the entire background geometry is 
important. In order to be able to apply the usual prescription of AdS/CFT to 
compute correlation functions, etc. it is desirable to construct a solution that 



■^*We thank Yaron Oz for pointing out the relevance of this work and fruitful discussions on this 
issue. 

■^^This approximation is viable only if £s/£ is small so that neglecting the higher derivative terms 
can be an approximation. As we mentioned several times this is not a tunable parameter in our 
set-up but rather should be determined phenomenologically by matching a corresponding quantity 
in the dual field theory. For example in certain models of holographic QCD p^ it turns out that 
this parameter is indeed small, about 1/7. 
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has the desired IR asymptotics ( p.22| ) and an asymptotically AdS UV geometry 
at the same time. This is done in section 6. One particularly interesting future 
direction is to find analogous solutions in JIB. 

The transition in the case Gd 7^ 

One can ask what happens if we turn on (small) gravitational interactions but 
ignore the higher derivative corrections. In this case the entropy of the BH 
solution (normalized by l/N"^) is given by, 

g(d-l)A($fe) 

where 6S is the correction that stems from the determinant around the saddle 
which is 0{Gd) with this normalization. The thermal radiation has a sin- 
gle contribution of the same order Stg = 0{Gd)- Therefore, the continuous 
transition would occur at $£ that is determined by 

g(d-l)A(<J.e) 

4GnN^ = Stg - SSsh- (8.2) 

There are two qualitatively distinct solutions to this equation. In the first case 
both sides of this equation may vanish. This is the same situation that we 
studied in this paper in the classical limit. If this is what happens then we 
conclude that the Gd corrections do not help avoiding the singularity and the 
transition still happens at the singular point where the scale factor vanishes. 
However, there is a second possibility, that the RHS is finite, hence the transi- 
tion is moved away from the singular point by the Gd corrections. In order to 
find which possibility is realized one should compute the spectrum of fiuctua- 
tions around the TG and the BH saddles at a temperature determined by $ft 
and see if the difference in the entropy can be matched to the LHS in ( ^.1| ) at 
a finite $£• One should also check that this $c does not lie beyond the point 
where the Jeans instability destroys the TG saddle. 

Adding other fields 

One can extend the analysis by including other fields in the action ( p.l|) . A 
particularly interesting case to consider is multiple scalar fields. One can ask 
what conditions should be imposed on the scalar potential l^($i, $2, ■ ' ' ) ^^ 
order to obtain continuous Hawking-Page transitions. 
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APPENDIX 

A. The scalar variables and reduction of the Einstein-scalar 
system 

The scalar variables of the Einstein-scalar system is defined as, In case of a 
neutral black-hole in d + 1 dimensions that is solution to ( p.ip , one only needs two 
such variables, that we call X and Y: 

where the function g is defined as g = log / and the constant 7 is given by. 



7^^/^. (A.2) 



These functions satisfy: 



§^-0^(1 -A- + 7)1 (A,4) 

This second order system is sufficient to determine all of the thermodynamic proper- 
ties (and dissipation) of the gravitational theory P . This is a reduction of the fifth 
order Einstein-scalar system to an equivalent second order system. 

It is straightforward to show that these equations combined with the following 
three. 



dA 
dr 




(A.5) 


d^ 

dr 


_le^o^X(<l>)e"^^*o(^W-^^"W)'^*, 


(A.6) 


dg _ 
dr 




(A.7) 



solve the original Einstein equations. Here $0 is a cut-off near the boundary. One 
can also express g and A in terms of the scalar variables directly from the definitions 
Ol): 

A{^)=A, + ^J idl>, (A.8) 

9^=1 f ^d^- (A.9) 
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Another useful equation relates the scalar potential to the scalar variables, that 
follows from ( [A.3|) : 



F($) = ^^^ (1 + r - X') e-2^/-\.(^W-^)'^*. (A.IO) 



The precise form of the overall coefficient follows from inserting ( |A.5[ ), ( |A.6| ) and 
( [A.7| ) in the Einstein equation ( |2.15| ). 



The zero T solution (thermal gas) corresponds to the Y = fixed point of eqs. 
O and (Q: 



^ = -7(l-X^)(l + -^^^^). (A.ll) 

The solution to this equation determines the zero T scalar variable Xq. In ||^ it is 
shown that, 

lim X{^) =Xo(<l>). (A.12) 



Equation ( |A.11|) provides a different representation of the scalar potential: 



V{^) = ^^^ (1 - X2) e-^^^-*o^oW)'^*. (A.13) 

B. Derivation of thermodynamics from the scalar variables 

The thermodynamics of the black-hole can directly be determined from the solution 
in terms of the scalar variables. The free energy is given by 

F($,) = ^ r d^h e('^-i)^(*'^) ^. (B.l) 



4^Gd i-j., d^h 

Our backgrounds satisfy the 1st law of thermodynamics S = ~dF/dT. Equation 
( p.l7|) directly follows from integrating this equation by noting the equation (|2.11| ) 

S = 47re('^-^)^("'^) + C(X-2), (B.2) 

and ignoring the l/N"^ corrections. The temperature in the scalar variable systems 
is obtained from 

T{^h) = ^^J_ ^^ e^(*'') V{^,) e^^^o ^(*) '\ (B.3) 

Derivation of this equation is more involved. It follows from the generalization to 
arbitrary dimensions, of the "useful relation" equation (7.38) of |^. The computa- 
tion is analogous to the appendix H.ll of [0 with slight changes due to different 
multiplicative factors in the equations in app. (|^) above. One finds the generalized 
"useful relation" as: 

S 1 /47r(rf-l)\'-'e->('^-i)/*^^ , , 

' ' (B.4) 



from which (|B.3|) follows immediately by use of 
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C. Sub-leading corrections to the geometry near T^ 

The sub-leading corrections are found by expanding the scalar variables as 

x{^) = -i= + 5X{^), Y{^) = Fo($) + 5r(<i>), (c.i; 



and inserting in (|3.9|) and (|3.10|). Then one obtains a linear differential equation 



for 6X that is solved by fixing the integration constant with the requirement of no 
singularity at $/»: One finds, 



Ck, -K$l-e-(''+'')(*h-*) 



5X($) - <( g^l^ L^(%J_4( / J-„-i ,.' (C.2) 

where we defined, 

7(^-1) ^ n d-^ /poN 

and the coefficients C, k and a are defined in (|3.1|) and ( p.2|) . The expression for Y 
follows from the leading term in ( p.lO ): 



yo($) = "^ ^ -^ . (C.4) 

One can also obtain the differential equation for 5Y from the sub-leading term in 
( p.lO|) but we shall not need it here. The correction term in ( p.l4| ) follows from the 



expression for 5X above. One obtains, 

_^ -K$ _)_... case i 

'-+7 ' 

C 
2 

in the limit $/» — ?■ oo. 



^ c$-° + ... case ii. 



D. Expressing Tc and $oo in terms of V^ 

Using equation ( p.3|) , we express As as. 



00 



.2A.(*,) _ ^2A($,) + ^*, ^ ^^ (Md--r)_\ ^_£_$,_2^J*.x_ p_^^ 



We consider the large $/j limit of this equation. The potential V^($) in this limit, is 
given by ( p.22| ). On the other hand, we can represent the potential V^($) (for any 
$) by equation ( [A.13| ). In the limit $/i — > oo the function X asymptotes to — l/v" 
as in ( |C1| ). Thus, by comparing ( |A.13| ) and ( p^.22D in this limit one finds^^. 



1)2 I 5 



e-^-^o'^^ = V;,_A_eV^*^ (D.2) 
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Using also the fact that the functions X($) and Xo($) become the same in the limit $/j -^ oo. 
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In the NCST case the normahzation of $ is given by ^ = A/{d—l). Then, substitution 
of (p.2| ) in ( p.l|) , use of ( p.22|) and the fact that T —i- Tc in this hniit one finds, 



On the other hand we know that As{^h) — !■ in this hmit. Thus, one finds. 



(D.3) 



Te 



^ 



4:71 



or 



V^o = (47rT,)2. 



(D.4) 



(D.5) 



To express $00 in terms of VJ30 we use the e.o.m. for $, ( |A.6D . We first recall 
that the constant $00 is defined as. 



lim $(r) = $00 f- 



(D.6) 



Use of (Q in (|70|) gives 



dr £7 



(D.7) 



Now, using ( p.2|) , (|C.1|) and the fact that ^ = 4/((i — l),in the NCST case, one finds 
in the limit $/i — ;■ cxd. 



lim — — 



$. 



7314 



$=* 



h 



7(rf-l) 
Using the definition ([A.2| ) and ,^ = 4/((i — 1), one finally arrives at. 



^^ r 



(D.8) 



(D.9) 



From this one can also express A^o in terms of Voo- We recall the definition of 



A. 



lim A{r) = —Aoo r. 

r— >oo 

Comparison of ( |2.21| ) and ( p.lO| ) gives. 



Aoo = ^^oo = ^^ = ^^ 
d — 1 d — 1 d — 1 



T 

-'- CI 



where the last equation follows from (U.5) 



(D.IO) 



(D.ll) 
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E. Curvature invariants 



Einstein's equations that follow from ( p.l|) read 



R,u - \g,uR = e (d,^d>'^ - ^g^Ad^A + ^g,uV{^). 



Tracing both sides one obtains, 



d+l 
d-1 



vm. 



(E.l) 



(E.2) 



One can express the right hand side in terms of the scalar variables of Appendix ^ 
Substituting ( [OD and ( |AlO|) in O one finds, 



i?-^ 
^-P 



{d - 1)X^ -{d+l){l + Y - X'^)e''-^o X 



-27/0* X 



(E.3) 



Both for the TG solution and for the BH solution in the limit ^h ^ 1 (which 
corresponds to the phase transition region T ^ T^) one has, Y ^ 0, X ^ Xq for all 
$ ^'^ . This is indeed the limit that we are interested in. Equation ([E.3|) becomes. 



d 



R^-[2dX',-id+l)]e 



P 



-2^foXo 



(E.4) 



The approximation becomes better and better as $/i gets closer to infinity in the BH 
case. In the TG case ( |E.4|) is exact. We are interested in how R scales with $ for 
large $. The biggest contribution come from the large $ region where we can use 
the asymptotic expression for the function X, ( |[J.1| ): Xq — )■ —l/y/d. Substituting 
this in ( [E.4|) and using ( p.2|) (X($) and Xq{^) become the same functions in the 
limit $/i S> 1), one finally arrives at. 



R 



dVr. 



d-l 






as $ > 1 



(E.5) 



To repeat, in the BH case the approximation is valid in the region $/i ^ 1. We see 
that R is proportional to V in the IR region. The same is true for the dilaton invariant 
[d^y. One can also show that the higher order invariants such as R^ua/sR'^""'^ also 
scale as (powers) of V in the IR. 

One can calculate the string frame invariants from the Einstein-frame ones by 
passing to the string frame metric. 



9s,^ 



jjtv 






Ricci scalar transforms as, 
R, = e ''-I 



R + 4^9^9^$ + A^^d^^d^^ 



(E.6) 



(E.7) 
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Except at strict $ == $/j where in the BH case Y has a singularity at this point. 
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where the indices are raised and lowered by the Einstein- frame metric on the RHS. 
Let us specify to the case of NCST normahzation ^ = 4/((i — 1) where also the 
relations ( p.28| ) hold. Using ([E.5| ) and ( p.28p we see that the first and the last terms 



in the square brackets cancel in the leading order, wheras the second term scales 
exactly as the subleading pieces in the first and the last terms. Therefore one finds, 

Rs oc d^6<l>{r), (E.8) 

where 5$ is the subleading piece in the IR: $ -^ $00^ + S^{r). In order to find this, 
one works out equation ( p.7|) to sub-leading order by substituting X in (|C.1|) and 



using ( p.2|) in the limit •I'/j — )■ 00. One then finds, 

(5$ oc 6^'^*°°'^ (E.9) 

in case i and, 

(5<l>ocr-" (E.IO) 

in case ii. Substituion in ( |E.8| ) finally yields, 

/2s oc e""*, {case i); $"°, {case ii), (E.ll) 

in the limit $ — )■ 00. We see that the string frame Ricci scalar vanishes near the 
transition region. The same is true for the dilaton invariant (5$)^ in the string 
frame. Higher order invariants vanish with a stronger rate. This shows that the a' 
corrections in the string frame become unimportant near the transition region. 

F. Bulk fluctuations and dissipation 



In this section, we present the details of the calculations in section p^. We first 
derive the fluctuation equation (|5.8|) that we reproduce here: 



(F.l) 
where prime denotes the derivative w.r.t. $ and the metric functions are defined in 
the <l>-frame: 

ds' = e^^(*) {-m)dt^ + dxU) + e'^^'^^y (F.2) 

One fluctuates the Einstein's e.o.m. with the metric fluctuations: 



g,, = e'''{l + -he'''^') , (F.3) 






g<S,<!, = —r- ( 1 + 2^**^ 
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and (5$ = 0. The last choice corresponds to the gauge (|F.2|) . We need only two of 
the zeroth order Einstein eqs. in this gauge: 

/" + dA'f - B'f = (F.4) 

A" - A'B' + -p— = 0. (F.5) 

The fluctuation equations for /iqo, h^ and /i$<j, are obtained by expanding the Einstein 
eqs. Ci = T^ to order e where the stress tensor at this order reads: 



rpv 



-ie-'^fh^^e-^-' ^5151 - ^5;:) . (F.6) 



One also needs to expand the Einstein tensors to this order. This is done with help 
of Mathematica. We shall present the final results. From the first order terms in e 
in the Einstein's equations one finds the following e.o.m for the fluctuations: From 
the {$, $} Einstein equation one finds, 

From the {0, $} Einstein equation one finds, 

h^^ = -^ (^K - (jh^ . (F.8) 

Finally, from the {i,i} Einstein equation one finds a complicated expression that 
involves hb, h[ and h'l that simplifies upon use of ( |F.7| ), ( [!''. 8| ) and the zeroth order 
eqs. ( |1''.4D and (^^ and reduces to ( |1:''.1| ). 



Now, we move on to the details that lead to ( p.7|) that we reproduce here: 



ImJ= ^e^e^^-^ [kK - hK') . (F.9) 



This is obtained by expanding the Lagrangian ( p.l|) to second order in e in fluctua- 
tions, as explained in [|19|, as follows: 

Define the fluctuation vectors h = (/loo, ^fe, ^<j>$)- Then the C(e^) term in the 
Lagrangian is of the form, 

C = -h^ Mu h + -h'^ M$$ h' + h'^ A/U h + -h^ M h+--- (F.IO) 

where the ellipsis refer to total derivative terms in $ and in t. Now, one formally 
promotes h to complex valued fields and writes [h should be understood as a complex 
vector in the following): 

C = h^' M^^h' + h^' M^h + --- (F.ll) 
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where remaining terms are unimportant as they do not contribute to Im J which is 
the conserved flux under the U{1) rotations h — )■ e'^/i: 

£ = 5$J+---, J = h^A'U^h' + {M^ + ---) h). (F.12) 

Here the elhpsis in C refer to terms that vanish on shell and the ellipsis in J refer to 
a symmetric matrix that vanish in the imaginary part of J hence irrelevant for our 
purpose. 

Thus the entire calculation is reduced to the computation of the 3x3 matrices 
M<j,$ and M<j,. Clearly, only the ^/gR piece of ( |2.1| ) will give rise to terms with 
single or double derivatives in $, as in (|F.10|) . Moreover, we are not interested in 



the t-dependence, therefore one can ignore the time dependence in the fluctuations 
( [F.3D . One expands ^/gR with help of a symbolic computer code and obtains M$$ 
and M$ as (after integration by parts in the double derivative terms and upon use 
of the zeroth order eqs. ( |F.4|) and ( [F.5|) ): 



dA-B I 1 d-2 n \ S/T l*^ 1) „dA-B 




M<j,$ = (d-l) f e^^-^ U ^ , M^ = ^—^ e' 



-2fA' * 2fA' 

[r^^f + 2dA'f * ; 



where the order of the components are according to h = (/iqo, hb, /i##) and the stars 
refer to terms that do not contribute to the imaginary part of the flux. One calculates 
the flux from (|1''.12|) . It reduces to (|1'\9|) upon use of the first order equations ( lb'. 7] ) 
and 
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